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ABSTRACT 

The sense of line asymmetry of solar p-modes in the intensity power spectra 
is observed to be opposite of that seen in the velocity power spectra. Theoretical 
calculations provide a good understanding and £t to the observed velocity power 
spectra whereas the reverse sense of asymmetry in the intensity power spectrum 
has been poorly understood. We show that when turbulent eddies arrive at the 
top of the convection zone they give rise to an observable intensity fluctuation 
which is correlated with the oscillation they generate, thereby affecting the shape 
of the line in the p-mode power spectra and reversing the sense of asymmetry 
(this point was recognized by Nigam et al. and Roxburgh & Vorontsov). The 
addition of the correlated noise displaces the frequencies of peaks in the power 
spectrum. Depending on the amplitude of the noise source the shift in the 
position of the peak can be substantially larger than the frequency shift in the 
velocity power spectra. In neither case are the peak frequencies precisely equal 
to the eigenfrequencies of p-modes. We suggest two observations which can 
provide a test of the model discussed here. 


Subject headings: Sun: oscillations; convection; turbulence 
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1. Introduction 

Duvall et al. (1993) found that solar p-mode line profiles are not exactly Lorentzian in 
shape but are somewhat asymmetric about the peak of the power. They showed that the 
power in their observed surface velocity spectra fell off more slowly at the low frequency 
side of the peak compared to the high frequency side of the peak. They also claimed that 
the sense of asymmetry is opposite in the intensity power spectra i.e. there is more power 
at the high freqnency side of the peak compared to the low freqnency side. The results of 
Duvall et al. have been recently conhrmed by a number of independent observations. 

Theoretical calculations have little difficulty in reproducing the observed velocity power 
spectrnm. According to the simplest models one expects the shape of velocity and intensity 
power spectra to be almost identical, as long as both of these are observed at the same 
place in the solar atmosphere, since the velocity and temperature eigenfunctions of a mode 
are linearly related to each other. Rast & Bogdan (1998) point out that nonadiabatic effects 
are unlikely to modify this relation. Nor does it appear that the process of the formation 
of absorption lines in the optical spectrum could cause the asymmetry of intensity power 
spectrum to be reversed (Sasselov, personal commnnication). It is therefore puzzling that 
the intensity power spectra are observed to be very different from the velocity spectra. 
Nigam et al. (1998) snggested that correlated noise added to p-mode oscillation amplitude 
affects the shape of the observed spectrum and could be responsible for the difference 
between velocity and intensity power spectra; they snggested that granules could provide 
the correlation. Roxburgh & Vorontsov (1997) suggested that velocity variation associated 
with acoustic sources, i.e., the turbulent eddies, are observed when eddies arrive at the top 
of the convection zone and this variation should be correlated with p-mode oscillations. 
Roxburgh & Vorontsov chose, arbitrarily, the relative amplitnde of the velocity variation 
associated with oscillation and with the noise sonrce to be of order unity and considered 
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them to be almost perfectly correlated in order to produce the reversal of line asymmetry. 

We calculate the relative contribution of turbulent eddies to the observed flux variation 
and its correlation with p-mode oscillation in a self consistent manner and show that this 
process explains the observed reversal of asymmetry in the intensity power spectra. The 
calculations are presented in section 2. The conclusion and some predictions are contained 
in section 3. 


2. Intensity oscillation and line asymmetry 

We hrst consider a simple one dimensional square well potential model, same as that 
considered in Abrams & Kumar 1996, to calculate the power spectrum of p-mode oscillations 
and include the contribution of fluctuating sources, which excite these oscillations, to the 
observed flux variation amplitude: 


dr'^ 


+ i2uj'y'ip^+ u^-V(r) 6{r - r^), 


( 1 ) 


where is the temperature perturbation function, is the quadrupole source amplitude, 
7 is the damping constant which is taken to be independent of r, and V (r) is a simplihed 
effective potential constructed to describe the solar p-modes: V = oo for r < 0, for 
0 < r < a the potential is zero, and V{r) = a constant, for r > a; a is the sound travel 
time from the lower to the upper turning point of a given mode and a is the acoustic cutoff 
frequency at the temperature minimum. The observed flux variation is roughly proportional 
to the Lagrangian temperature perturbation at the place in the atmosphere where lines are 
formed|^. 


perturbation of the convective flux by oscillations could however modify the relation 


between the observed flux variation and the Lagrangian temperature perturbation. 
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For sources lying in the evanescent region i.e. r* > a, the amplitude of the temperature 
variation observed at point Tq can be easily calculated by solving Eq. (1) and is given by 


Tuj{ro) = S^G^{rs, Vo) + exp(i 0 - 77 ), 

1^2 


( 2 ) 


where 


Gco(rs,ro) = 


exp[—cc; 2 (r’ — a) 


UJ2 


ijj2 sin(a;ia) cosh(a;25r) + cui cos(a;ia) sinh(cc;25r) 


(3) 


102 sincuio + uji coscuia 

is the Green’s function for Eq. (1), ui = + * 270 ;)^/^, U 2 = — u? — 

6r = Ts — a, and 7 is the optical depth between Tq, the point where the observation is made, 
and the top of the convection zone. The dimensionless parameter f3 is the ratio of the 
observed amplitude of the flux variation of an eddy to that of the oscillation, and 0 is its 
phase difference with respect to the mode it excites at radius Tq. 


The real and the imaginary parts of the Green’s function have very different shapes as 
shown in Eg. 1, and this enables a correlated noise source to modify the observed shape of 
the power spectrum as discussed below. 


In the absence of the contribution of the noise source to the observed flux variation, 
i.e. when 0 = 0, the shape of the power spectrum is shown in the left panel of Fig. 2. This 
case should correspond to the velocity power spectrum of solar p-modes since the random 
component of the velocity held in the solar atmosphere — where optical observations are 
made — is believed to be uncorrelated with the turbulent eddies in the convection zone 
which excite p-mode oscillations. The asymmetry in our calculated spectrum does in fact 
look very much like that seen in the velocity power spectra of solar p-modes. 


Unlike the velocity variations, the hux variation associated with turbulent eddies are 
observed when the eddies arrive at the top of the convection zone. These variations are 
correlated with oscillations these eddies generate; this correlation is parameterized above 
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by (3 and 0. The addition of this ‘noise’ term to the oscillatory signal modifies the relative 
contribution of the real and the imaginary parts of the Green’s function to the observed 
intensity spectrum (which is proportional to |Tt^(ro)p) and thereby modihes the shape of 
the p-mode lines in the intensity observations (see Fig. 2, right panel). We hnd the reversal 
of asymmetry over a wide range of (3 and 0 as shown in £g. 3; these line profiles are 
similar to the observed reversal of asymmetry seen in the solar intensity data. Note that 
the background in the intensity spectrum is significantly non-zero, and in fact if this model 
is correct we expect the reversal of asymmetry to disappear when the background in the 
spectrum vanishes. 

The addition of the correlated noise source to the observed flux variation significantly 
modifies frequencies of peaks in the power spectrum as is clearly seen in Fig. 3. A change in 
the amplitude of the noise source, as quantified by the dimensionless parameter (3, from 0.5 
to 2 shifts the peak frequency of 2.4mHz mode by about 0.4 /iHz. The peak frequencies in 
the intensity power spectra are different from the peaks in the velocity spectra and neither 
of these are the true eigenfrequencies of p-modes. The frequency shift in the intensity 
spectra depends on the amplitude of the background in addition to the effects that cause 
shift in the velocity spectra (cf. Abrams & Kumar 1996, Christensen-Dalsgaard et al. 
1998), and so it is safer to use the velocity spectra to determine the true eigenfrequencies of 
solar p-modes. 

We calculate the parameters (3 and 0 for a realistic model below. 


2.1. Relation between wave and source amplitudes 

We consider a plane parallel atmosphere which sits in a constant gravitational held, 
g, and consists of two semi-infinite layers, the lower layer being adiabatic and polytropic. 
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and the upper isothermal. The pressure, p, density, p, and temperature, T, are continuous 
across the interface between the two layers. In the lower layer the adiabatic and polytropic 
indices are related by 7 = 1 + 1/m. The z coordinate measures depth below the level at 
which the adiabatic layer would terminate in the absence of the isothermal layer. In the 
adiabatic layer the thermodynamic variables have a power law dependence on the depth 
z e.g. p = pt{z/, p = pt{z/zt)^, the sound speed = gzim, and the scale height 
H = z/{m + 1), where the subscript t denotes quantities evaluated at the top of the 
adiabatic layer. In the isothermal atmosphere p and p are proportional to exp(^/ifi). 

Following Goldreich and Kumar (1990), hereafter referred to as GK, we use Q = Pi/p 
as the wave variable, where Pi is the Eulerian pressure perturbation. The inhomogeneous 
wave equation in the adiabatic atmosphere is 

^ ^ 1 ^ _ 2 

dz"^ dz (? dP ^ dP\Cp) pdz\Cpj 

where s is the Eulerian entropy perturbation associated with the turbulent convection. 
The source term S consists of a monopole and a dipole piece which arise as a result of 
volume change due to fluctuation in entropy and buoyancy force variation associated with 
this volume change respectively. For simplicity we have not included the Reynold’s stress 
source term associated with fluctuating internal stresses in the convection zone. Inclusion 
of Reynold’s stress makes the algebra considerably more tedious but does not modify the 
main results of this section in any essential way .0 

We expand the enthalpy wave function Q in terms of the normal modes, Qa{z), which 
are normalized to unit energy, 

^ GK showed that the monopole and the dipole source terms in S cancel leaving the 
remainder which has magnitude of same order as the quadrupole Reynold’s term. 
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[AaQa exp{-iut + ikh ■ x) + A* Q* exp{iut - ikh ■ x)], 


( 5 ) 


where A is the horizontal cross section area of the atmosphere, and kh is the horizontal 
wave vector. Substituting this expansion into equation (4) and multiplying both sides with 
Q’^exp{iujt — ikh ■ x), and integrating over space and time, we obtain 


Aa{t) = Sexp{iut - ikh • x). (6) 

Using the wave equation (4) this can be rewritten as 

The change to mode amplitude, SAa, due to a single eddy of size h over its lifetime t/j, 
which is located at xq and lived at time can be easily estimated from the above equation 
and is given by 


r/i u\ PC^h^^hSh{to)d^Ql ^ u ^ \ : -^1 

dAAt) ^ - 2icc;c ^U2 -exp[*ct;„to - ■ xo - 7«U “ ^o) + i<Pi\, (8) 

where the phase 4>i depends on the temporal and spacial properties of the eddy, and the 

factor exp[— 7 Q,(t — to)] has been multiplied to model wave damping that causes the change 

to mode amplitude to decrease with time for t > GK showed that the resonant scale 

height size eddies are most efficient at exciting low frequency modes (modes of frequency 

less than about 3mHz); for these eddies h ~ if and th ~ The depth {z^j) of these 

eddies below the photosphere, or the isothermal atmosphere in our model, can be estimated 

from the constancy of the convective flux pv^ 


MtUJac 


3/(m+3) 


■} 


Zca ~ Zt 


00 


( 9 ) 



9 


where Mt is the Mach number of the turbulent convection at the top of the convection 
zone, and ujac is the acoustic cut off frequency in the isothermal atmosphere. 

The change to the enthalpy perturbation associated with a mode a, SQa(zo), in the 
isothermal atmosphere at location Zo due to the eddy is given by 


^Qai^Zoj t) 


6 AaQa{zo) exp{-iuat + fkh ■ x) 


Substituting for from Eq. (8) we hnd 


( 10 ) 


exp 


T(z^ 


dz"^ 


■(7a+ici;«)(f-to)+i0i+*kh ■ (x - xq' 

(11) 


where ST^d = Sh/cp is the temperature fluctuation of the eddy. Low frequency p-modes 
are most likely excited near the top of the convection zone in a region where waves are 
evanescent (cf. GK). The wave function is nearly constant in the evanescent region and 
its derivative d’^Qajdz^ ~ {uj"^/g‘^)Qa- Moreover, the normalized eigenfunction in the 
evanescent region is (see GK) 


Qc 




p]/‘^g{ra-2)l2 

Making use of these equations, the expression for 5Qa reduces to 


( 12 ) 


f) 


Ag^ 


Cp5T^d{z^Ao) - to) + *kh ■ (x - Xq) + i0i]. (13) 


Finally substituting for z^^ from Eq. (9) and using the relation gjzt zz the above 
equation further simplihes to 
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'tkh { uj \ ^ 

5Qa {zo, t) - t^MI — Cp6Ted{zuj, to) exp[-( 7 a + iua) (t - to) + ikh ■ (x - xq) + i0i]. 

vA \ OJac / 

(14) 

We next consider the relation between 6Q and the associated change in the temperature 
amplitude of the mode at the height in the atmosphere where optical observations are made. 
For the nearly isothermal solar atmosphere the Lagrangian temperature perturbation and 
the pressure perturbations are related by 


AT = Ap 


© 


Ap 

pCp 


Pi + pgZ 

pCp 


(15) 


where is the radial displacement associated with the wave, and for low frequency acoustic 
waves it is given by (see GK for detail) 


1 

dQ {i-l)Q 

Q 

1 

1 

e 

to 

, ,2 
^ac 

dz '-)Hi 


. 7 40^acJ 


Therefore, 


(16) 


AT 


Co- Q 

^ac 


(17) 


Substituting this into Eq. (14) we obtain the change to the Lagrangian temperature 
perturbation in the isothermal atmosphere as a result of wave generated by the entropy 
fluctuation of an eddy 


ATa{zo,Xh,t) ^ iF{u, kh)STed{z^,to) exp[-{'ya + i(^a){t - to) + fkh ■ (x - xq) +# 1 ], (18) 


where F{uj, kh) is a dimensionless function dehned by 
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F{u, kh) 


khZt 

A 



Ml 


(19) 


The total observed temperature fluctuation in the atmosphere, as seen in a line or 
continuum, at some location on the solar disk, consists of the contribution from oscillations 
and the attenuated emission from turbulent eddies at the top of the convection zone. 


5Ta{zo,yihA) = i -F(i^,/i^/i)^T’erf(2:^,Xo,fo)exp[-(7„+ icn«)(t-to)+*kh • (x-xo)+i0i] 

eddies 

+(5Ted(2:t,x,t)exp(-?7), (20) 

where STed{zt,^,t) is the temperature variation associated with eddies at the top of the 
convection zone, and rj is the optical depth of the line, used in the observation, at the base 
of the photosphere. 

The summation above is carried out to include all resonant eddies. To project out 
modes of a particular horizontal wave vector we multiply both sides by exp(—ikh ■ x) and 
integrate over the area A of the box. Individual modes are isolated by Fourier transforming 
in time and the result is 


kh, cn) 


-A 

{uj - UJa) + il 


- Y. F{u,kh)STU 

^ eddies 


Xo, to) exp{iujto - fkh ■ xq + # 1 ) 


+FZt Y ^Ted{zt, Xi, ti) exp{iuti - ikh ■ xi + - 77 ), (21) 

eddies 


where the phase depends on the details of the spacial and temporal structure of eddies, 
and it should be roughly equal to the phase so long as the temporal property of the eddy 
does not change dramatically between the depth where the waves are excited Y) and the 
top of the convection zone. 
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The observed power spectrum P{uj) = |(5Tap. Assuming that ST^d for different eddies 
are uncorrelated, we retain only those cross terms in which correspond to the same 

eddy; (to — ti) in this case is equal to (z^j — Zt)/ved = SUt is the time for downward moving 
eddies to travel from the top of the convection zone to the depth where they most efficiently 
excite the mode. The power spectrum at the frequency Ua is thus given by 






P(ua) ~ 2 3 kh)STedj^ - ^F(ua, kh) sin(uAttt) \5TedY exp(-? 7 ), (22) 

Ztlin 7a 

In deriving the above equation we have assumed that the duration of the observation is 
greater than the mode lifetime; for smaller observational time T, the second term is reduced 
by a factor of 


Due to the multiplicative factor sm(uJaStu) the second term in equation (16) averages 
to zero if downward and upward moving eddies contribute equally to wave excitation. 
Thus the reversal of asymmetry seems to require that downward moving eddies contribute 
more to the generation of acoustic waves; some recent observations suggest that modes are 
probably excited by the downflows (cf. Goode et ah 1998). 

Comparing the two terms in Eq. (21) with the corresponding terms in equation (2) 
we find that the phase difference between the observed flux variation of a mode and the 
turbulent eddy is 0 = tt — uStu] In deriving this result we expanded the Green’s function in 
equation (3) in the neighborhood of an eigenmode: ~ ui[((jJ — uia) + i 7 a]~^/( 8 aQ;^). The 

value of 5ttt is shown in fig. 4, and we see that 0 ~ 1 rad for low frequency modes. 


Moreover, from equations (21) and (2) we find 


/9 


^oTt 


ui„ 


2m—1 


8(auJac)khZt \uJo 

At the top of the solar convection zone A4t ~ 0.3, Zt 


M 


-3 


(23) 


300km, g ~ 2.77x10^ cm s 


-1 
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m ~ 1.5, Tt ~ 100s, and in the solar photosphere Uac ~ 3.3xl0“^rad s“h The sound 
travel time from the lower to the upper turning point, a, for this model atmosphere is 
approximately uja/{gkh), and therefore akh ~ uja/g- Substituting all this into equation 
(23) we hnd the value of f3 to be about 100 for a mode of frequency 2mHz, and this is 
independent of the horizontal wave number kh- A more realistic model of the sun considered 
in the next subsection yields essentially the same value for (3 when sources are taken to 
he at a depth ~ 100km as suggested by the mixing length theory (see eq. 9). For 
this source depth the line asymmetry in intensity power spectra is reversed as long as the 
effective optical depth, at the top of the convection zone, for the wavelength band used in 
the intensity measurement is less than about 5 so that f3' = /5exp(— 77 ) is greater than ~ 0.5 
(see £g. 3). 


2.2. Results for a solar model 

The equation for (3 for a solar model, corresponding to equation (23) of the plane 
parallel atmosphere, can be derived in a manner analogous to that given in § 2.1 and so we 
outline a few of the main steps in the derivation. 

It is convenient to work with the displacement wave function in this more general case, 
instead of the enthalpy perturbation used in § 2 . 1 , which we expand in terms of the complete 
set of eigenfunctions ^cX^m 


I 


1 

71 


^ AaicXtra exp(-fct;«f) + C.C. 
a. 


(24) 


The change to the radial displacement amplitude of a mode, due to one eddy, at radius Tq 
in the solar atmosphere is given by (see Goldreich et ah 1993) 



14 


xc ^ ^ ~ (9p\ ^ 


'8Q 


dr 


Y, 


im' 


(25) 


The notation used here is same as in equation (11). The change to the Lagrangian 
temperature perturbation in the atmosphere resulting from the excitation due to one eddy 
is calculated using the above equation and is given by 


AT„(r„ 




2(7 - l)7b 





dr 


YL^T^, 


(26) 


where Hi is the pressure scale height at the place in the atmosphere where the optical 


observation is made, is the scale length of resonant eddies, and is the radius where 
most of the wave excitation takes place. It is straightforward to calculate (3 from this 


equation 




2(7 - l)H^Ur^) 

[dp, 




'CQr(^o) 




dr 


-1 




aUJn 


(27) 


where the subscript t refers to quantity evaluated at the top of the convection zone, and a 
is the sound travel time from the lower to the upper turning point of the wave. Figure 4 
shows graphs for (3, for three p-modes of a solar model, as a function of depth below the 
photosphere. The value of (3 for low frequency modes, at reasonable excitation depth, is 
consistent with our estimate of § 2 . 1 . Also shown in £g. 4 is the eddy travel time, 5tu, 
from the top of the convection zone, from which we can estimate the phase difference 
(j) = n — UaStu- The values of (3 and (j) for the solar model seem to lie within the parameter 
range needed for the reversal of line asymmetry. 


Specihcally, if a mode of frequency ~ 2 mHz is excited by quadrupole sources lying 
about 100 km below the top of the convection zone then (3 ~ 100 (see £g. 4); this is same as 
estimated in §2.1. For flux variation observations using the central part of the Ni line, one 
of the hlter bands used by the GONG and the MDI instruments, the value of 7 , the optical 



15 


depth at the top of the convection zone, is about 6 . Thus jS' = /3exp(— 77 ) for this source 
depth and the optical line is about 0.3, and hgs. 3 and 4 suggest that the line asymmtery 
in the intensity power spectra should be opposite to the velocity power spectra. However, 
if the acoustic sources are located at a depth of 200 km or greater, as suggested by Kumar 
and Basu (1998), then (3 would be much smaller (see £g. 4) and the line asymmetry in 
the intensity power spectra obtained using the Ni line should be the same as the velocity 
spectra. The line asymmtery in intensity power spectra, for deeper lying sources, is expected 
to be reversed when observations are made in the continuum for which 77 ~ 1 . 

The asymmetry reversal requires that the effective observed amplitude of the noise 
source, f3' = /5exp(— 77 ), exceed some minimum value of order 0.1 (see fig. 3). Thus we 
expect the velocity and the intensity power spectra to have the same sense of asymmetry, 
irrespective of the source depth, in those observations carried out using optical lines formed 
high up in the chromosphere so that their optical depth 77 at the top of the convection zone 
is large. This prediction provides a verifiable test of the model discussed here. 

For a fixed source depth the value of (3 decreases rapidly with increasing mode 
frequency (see fig. 4). Thus if modes of different frequencies were to be excited at the same 
depth, the line asymmetry reversal seen in the intensity power spectra should disappear 
at some frequency of order 3.5 mHz. Asymmetry reversal in the observed spectra at this 
frequency suggests that the higher frequency modes are perhaps excited higher up in the 
convection zone. 


3. Discussion 

We find that the observed reversal of asymmetry in the intensity power spectra of solar 
p-modes can arise as a result of contribution to the observed flux from turbulent eddies that 
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excite these oscillations. The flux variations from these eddies is correlated with the p-mode 
oscillations they generate (as suggested by Roxburgh & Vorontsov 1997, Nigam et ah 1998), 
causing a change to the relative contributions of the real and the imaginary parts of the 
Green’s function to the power spectra and consequently to the observed line shape. We 
have calculated the correlation between the eddy flux variation and the oscillations and And 
it to be well within the range required to reproduce the reversal in the line asymmetry. The 
velocity power spectra of solar p-modes are unaffected by the contaminating noise source 
because the random velocity held in the solar atmosphere, where optical line observations 
are made, are completely uncorrelated with the oscillations. Intensity power spectra on the 
other hand contain some contribution from turbulent eddies at the top of the convection 
zone. These eddies excite solar p-modes as they travel downward, and the flux variation 
associated with them is correlated with the oscillation they generate. 

A reversal of asymmetry in the model problem considered by Roxburgh & Vorontsov 
(1997) was obtained when a superposition of dipole and quadrupole sources was considered. 
However, we And that asymmetry reversal can arise even when modes are excited by 
quadrupole or dipole sources alone. Furthermore, we have calculated the relative amplitude 
and the phase correlation in a self consistent manner, and And these to he in the parameter 
range that leads to a reversal of asymmetry. 

If the model considered here is correct then we expect the asymmetry reversal to 
disappear when the observed background in the power spectra is vanishingly small. 
Moreover, the use of spectral lines which are formed deeper in the solar atmosphere, to 
observe solar oscillations, should have larger amplitude for the flux variation associated with 
turbulent eddies at the top of the convection zone and hence the line asymmetry reversal 
should be more dramatic. Conversely observations made in optical lines formed high up 
in the chromosphere, so that their optical depth at the bottom of the photosphere is very 



large, should show intensity and velocity power spectra to have similar shapes. Recent work 
of Schou (1998) appears to support this prediction. 

We are grateful to John Bahcall for encouraging us to think about this problem, for 
several helpful discussions, and for his comments on the paper. We thank J. Christensen- 
Dalsgaard for snggesting several improvements to the paper and for mentioned to us the 
recent work of Jesper Schon. This work was snpported by a NASA grant NAG5-7395. SB 
was partially snpported by an AMIAS fellowship. 
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Fig. 1.— The real (left panel) and the imaginary (right panel) parts of the Green’s fnnction, 
for the simple 1-D model problem, at a hxed source and observer location as a function of 
wave frequency. The source is located in the evanescent region approximately 200 km above 
the upper turning point, and the observer is at a place corresponding to the temperature 
minimum in the solar atmosphere. 
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Fig. 2.— The left panel shows the velocity power spectrum for the simple 1-D model 
problem i.e. it does not include contribution from the noise source. The graph in the right 
panel includes the flux variation from turbulent eddies and should correspond to the intensity 
power spectrum; the values of the parameters /5exp(—?]) and 0 are 1.5 and 50° respectively. 
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Fig. 3.— Power spectra for several different values of the parameters (/3, 0). The continuous 
lines are line prohles when 0 = 0°, the dotted lines for 0 = 50° and dashed lines for 0 = 120°. 
For 270° > 0 > 90° the sense of asymmetry is same as in the velocity power spectrum (see 

fig- 2). 
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Fig. 4.— The left panel shows the value of (3 for three different modes of the sun as a 
function of depth measured from the top of the convection zone; the degree (£) of all the 
three modes is 50, and z/ = 2.08mHz for the solid curve, 2.67mHz for the dashed curve and 
3.2mHz for the dot-dash curve. (3 has a very weak dependence on it changes by about 
20% between £ of 2 and 50. The panel to the right shows the eddy travel time from the top 
of the convection zone {Stu) as a function of depth. 





